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Abstract. In this paper we use Lagrange- Poincare reduction to understand 
the coupling between a fluid and a set of Lagrangian particles that are supposed 
to simulate it. In particular, we reinterpret the work of Cendra et al. [CMR01 
by substituting velocity interpolation from particle velocities for their princi- 
pal connection. The consequence of writing evolution equations in terms of 
interpolation is two-fold. First, it gives estimates on the error incurred when 
interpolation is used to derive the evolution of the system. Second, this form 
of the equations of motion can inspire a family of particle and hybrid particle- 
spectral methods where the error analysis is "built-in". We also discuss the 
influence of other parameters attached to the particles, such as shape, orienta- 
tion, or higher-order deformations, and how they can help with conservation of 
momenta in the sense of Kelvin's circulation theorem. 

1. Introduction 

In this paper we seek to understand, from a geometric point of view, how a 
set of Lagrangian particles can be used as a computational device to numerically 
simulate an ideal fluid. Specifically, there are certain quantities which we can 
associate with a finite set of particles, such as their positions, velocities, or shape 
change to various orders (e.g., as an evolution in SL(3)). All these attributes may 
be derived by integrating the equations of motion. 

However, as we point out in this paper, much insight can be gleaned by ap- 
pealing to the process of Lagrange- Poincare reduction. Specifically, given an ideal, 
homogeneous, inviscid, incompressible fluid on M, the configuration space may be 
described by the group of volume preserving diffeomorphism, SDiff (M). If © is an 
iV-tuple of distinct points in M, then we define the isotropy subgroup 

G Q := e SDiff(M) I V(0) = ©}• 

The particle relabeling symmetry of the system allows us to project the equations 
of motion onto the quotient space T SDiff (M) / G & . Upon choosing an interpolation 
method (defined in §2.3p . we obtain an isomorphism to the vector bundle TX©g Q , 
where X is the configuration space of point particles and g© is a vector bundle 
over X whose fibers are isomorphic to the Lie algebra of the symmetry group. We 
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use the resulting Lagrange-Poincare equations to describe the coupling between a 
fluid and its computational particles in terms of interpolation methods; we also 
propose a new family of particle methods. 

1.1. Previous Work. It was shown in [Arn66j that the Euler equations of mo- 
tion for an ideal, homogeneous, inviscid, incompressible fluid on an oriented Rie- 
mannian manifold M are the spatial (or Eulerian) representation of the geodesic 
equations on the group of volume preserving diffeomorphisms, SDiff(M). This 
observation gave rise to a new perspective on fluid mechanics which lead to many 
developments, notably the proof of well posedness |EM70] and various extensions 
ranging all the way to charged fluids, magnetohydrodynamics, and even complex 
fluids with advected parameters (see, e.g., [Hol02j, [GBR09J). All of these system 
are Lagrangian on the tangent bundle of groups of diffeomorphisms of a Riemann- 
ian manifold M. Additionally, all of these theories utilize the particle relabeling 
symmetry of the system to perform Euler-Poincare reduction and bring the dy- 
namics to the Lie algebra of this group |MR99| chapter 13]. 

We may consider reducing by subgroups of the diffeomorphism group, and there 
do exist frameworks for accomplishing this. This would be a special case of 
Lagrange-Poincare reduction [CMR01J. In particular, we may consider reducing 
by isotropy groups of a set of points in M. Such an approach is already mentioned 
in |MD10j for the purpose of landmark matching problems; see also the references 
cited therein. However, to the best of our knowledge, Lagrange-Poincare reduction 
has not been performed on such systems in the framework of |CMR01j . 

1.2. Outline. In §2] we establish our notation and review the notion of a gener- 
alized connection (also called an Ehresmann connection [MMR90J) as described 
in |KSM99j . We will also state a few useful formulas which result from having a 
generalized connection. In §2.21 we will discuss diffeomorphism groups and quo- 
tients spaces. In §2.31 we will articulate the relationship between interpolation 
methods and generalized connections. The equivalence of interpolation methods 
and principal connections is given in Appendix |Aj Then, in §j3j we will describe 
the unreduced system for an inviscid fluid before reducing it by the isotropy sub- 
group of a finite set of particles. In §H we will discuss reduction by higher-order 
isotropics which will allow us to study particles with orientation, shape, and other 
attributes. Finally, in §51 we will formulate a family of particle methods induced 
by an interpolation method and discuss some implications for the error analysis of 
these methods. We will find that it is possible to construct hybrid particle-spectral 
methods for fluids within this family. Moreover, we will show that the vortex blob 
algorithm fits within this family of methods and that the horizontal equations are 
a guide for corrections that allow for the deformation of vortex blobs. We close 
with §|6l where we summarize how to extend these constructions to complex and 
other fluid models. 
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2. Preliminary Material 

Before introducing our contributions, we review generalized connections and 
volume-preserving diffeomorphisms and prove a few important theorems. 

2.1. Generalized Connections. In this section we introduce the notion of a gen- 
eralized connection, as presented in |KSM99j . and prove some useful propositions 
for the purpose of this paper. 

Definition 2.1. Let tce '■ E — » M be a vector bundle and te '■ TE — > E the tangent 
bundle of E. Then the vertical bundle is the vector bundle ttv(e) '■ V(E) — > E 
where V(E) := kernel(T7rs) and tiv(e) '■= t e\ V (e)- 



The vertical lift operator, v * : E © E — > V(E), is defined by 

(v m + ew m ) , 



v^{v m ,w m ) := 



e=0 

-1/ 



for all v m ,w m G E m := ix E (m). This establishes an isomorphism between the 
vector bundles proj x : E © E — >■ E and iiv{E) '■ V(E) — > E, where pro^ denotes 
the projection onto the first summand. The fiber derivative, j£ : E — > E*, of a 
function / G C°°(E) is defined by 



{fe (e) ' 6 ') := (dm^(e,e'))=f e 



UXe + ee')). 

=o 



We shall need a notion of partial differentiation with respect to the base space M. 
That is, we need to understand what is meant by J^j- evaluated at e G E, as an 
element of T*M over the base point m = tte{M). This is obtained by the use 
of a covariant derivative. For this, we define the vertical drop by v± : V(E) — > 
E := proj 2 o(f^) -1 , where proj 2 : E © E — > E is the projection onto the second 
component. 
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Definition 2.2. Let tie '■ E — >■ M be a vector bundle. A horizontal bundle is a 
subbundle, H(E) C TE, such that TE = H(E) © V(E). This defines projectors 
hor : TE — y H(E) and ver : TE — y V(E). The vertical projector, ver, is called 
a generalized connection. IfC(I;E) denotes the set of smooth curves in E on an 
open interval I C R, then the covariant derivative induced by H(E) (or ver) is the 
map : C(I; E) E given by 

De ( ( de 



—— = v\ ver . . 
Dt l \ \dt 

The choice of a generalized connection on a vector bundle tte '■ E — > M induces 
the horizontal lift operator : E@TM — y H(E) by the condition that h^(e, m) is 
the unique horizontal vector in the fiber T e E such that TiTE(h^(e, m)) = m, where 
e e E and m G T m M. This implies that 



(te®Tixe)\ h{e) 



Thus, the choice of a generalized connection introduces the partial derivative of 
the function / G C°°(E) with respect to m as a vector bundle map : E — > T*M 
over the base M defined by 

Qf \ 

(e),5m) := (df(e),h t (e,5m)) , 



dm 

for all e G E and 5m G T m M. 

By construction, this means that the total exterior derivative df acting on the 
velocity of a curve e(t) G E over m(t) G M be written as the sum 

de\ / 9f f . dm\ /df De 

" /(e) '^/ = W (e) '^/ + W (e) '^ 
If M is a Riemannian manifold we may use the Levi-Civita connection, de- 
fined on the vector bundle TM. The torsion free property of the Levi-Civita 
connection is equivalent to the statement that = T^, where are Christoffel 
symbols of the connection in an arbitrary local coordinate system. The covari- 
ant derivative with respect to the Levi-Civita connection is given (locally) by 

= it ^% vl ^& J '' wnere mj denote local coordinates of the point m G M. 
Therefore, the horizontal lift induced by the Levi-Civita connection is given lo- 
cally by 

(2) h\{m\ v j ), (m\ 5m ] )) = {{m\ v j ), (5m l , -rJ-uW)). 
This allows us to easily prove the following proposition. 

Proposition 2.1. Let a G Q l (M). A torsion free connection on TM (such as a 
Levi-Civita connection) induces the identity 

(3) da(m)(v,w) = /^(w),v\-/^(v),w 
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where v, w G T m M. 

Proof. This may be verified in a local coordinate chart, where a(m) = a.i(m)dm l . 
Viewing a as a function on TM we find that equation (J2J) implies 

Since da(m)(v, w) = (j^jj — v^w\ we can see that (J3J) follows from the torsion 

^k 

ij ji ■ 



free property, rj-L- = T k ji . □ 



This proof is extended to infinite dimensional manifolds by noting the irrelevance 
of the local coordinate description of rf-. It can also be extended to one-forms on 
M with values in a vector bundle E — > M by defining the exterior derivative on a 
tensor product by 

d(e <g> a) := e da 

where e is a section of a the vector bundle E — >■ M and a G f2 1 (M) is an ordinary 
one-form on M. 



2.2. Diffeomorphism Groups. Let M be a finite dimensional connected Rie- 
mannian manifold with metric (-,-) M '■ TM © TM — > R. The set of volume 
preserving diffeomorphisms ,SDiff(M), of M is an infinite dimensional Frechet Lie 
group. A tangent vector v v G SDiff(M) over a group element ip G SDiff(M) is 
a map : M — > TM such that v v (m) G T^( m )M. The Lie algebra jtdiv(-W) °f 
SDiff (M) is the set of divergence free vector fields on M. Therefore, T v SDiff (M) = 
{uo (p \ u e X div (M)}. 

Let © = (0i, . . . , ©at) G M w be an iV-tuple of distinct points in M. We define 
the isotropy group 

G Q := G SDiff(M) | V(0) = ©}, 

where SDiff (M) is acts on G by the diagonal action. It is elementary to 
see that the quotient of SDiff (M) by the right action of G Q is the set 

X := {(mi, . . . , m N ) G M N \ rrii ^ raj}. 

We shall define in £J3] a Lagrangian system on T SDiff (M), invariant under right 
translations by elements of G & and implement Lagrange-Poincare reduction along 
the lines of [CMR01J. This task is less trivial than it may first sound because 
(T SDiff (M))/G Q ^ T (SDiff (M)/G ). We denote elements of T SDiff (M)/G Q by 
[v v ] for each v v G TSDiff(M). We will find that (T SDiff (M))/G G is a vector 
bundle over X. Thus, Lagrange-Poincare reduction will provide dynamics which 
describe how the motion of a finite set of particles couples to the fluid. 
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2.3. Interpolation Methods as Generalized Connections. Define the vector 
bundle 

0© := {(z,0 I x e X, £ e q x } 

where q x is the Lie algebra of the isotropy group G x for each x G X, i.e., the set 
of divergence free vector fields which vanish at x G X C M . As a side note, 
the bundle $j is identical to the adjoint bundle SDlff ^ xg0 when equipped with 
the fiberwise Lie bracket [(x,£ x ), (x,r) x )] := (x, [£ x , ^Jjacobi-Lie) and the projection 
k(x,£ x ) = x. 

We will ultimately identify the quotient (T SDiff (M))/G with TX © g . 

Definition 2.3. A 3£di v (M) -valued one-form X : TX — > Xdi v {M) on X such that 
I(x)(x) = x for x G TX and x = Tx(x) is called an interpolation method. 



Figure 1. Schematic representation of an Interpolation Method 

Example 2.1. Initially, one would desire that j M {T(x)(m),^ x (m))dvol(m) = for 
all (x,£ x ) G g©. Such an interpolation method would correspond to choosing the 
"mechanical connection" (see [MMR90J section 2.4). However, this interpolation 
method does not exist. Indeed, if such an interpolation method did exist, then 
T(x) for some non-zero velocity x G TX would need to vanish everywhere, yet be 
non-zero at the locations of some of the particles. This would not define something 
which is smooth. However, if we choose an inner product which makes j£div(^) 
into a Hilbert space we may be able to do something more systematic. For example 
if M = R d , and I : X(R d ) -> X(R d ) is a positive definite SE(3)-invariant operator, 
we may define the inner product 





with associated reconstruction mapping 



X(x)(m) = ^G(| 



m — Xi\)x\ 
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where is Green's function of I on the space X(R a! ). For example, if I = 

1 — a- A for some a > 0, then = exp(—\x\/a). This is not directly applicable 

in the context of ideal incompressible fluids since this interpolation method does 
not produce divergence free vector fields and the operator 1 — a ■ A is not naturally 
identified with ideal fluids. However, this construction is natural in the case of 
the EPDiff equation with respect to the Lagrangian induced by the H 1 norm (see 
[HSS091 part II] and reference therein). 

Returning to the general case and using coordinates (x, £ x ,x, £») f° r Tg & , we 
can see that the choice of an interpolation method X induces a horizontal space 
for the vector bundle g given by elements of the form 

(x,£ x ,x, %,X(x)}) E Tg Q . 

In particular, we obtain the horizontal and vertical projections 

hov(x, £ x , x, 4) = (x,£ x ,x, %,X(x)}) 

ver(x, x, 4) = (^,^,0,4- [fx, 

Proposition 2.2. The covariant derivative induced by the horizontal space induced 
by an interpolation method, X, is given by 



D(x,Q ( ^ d£ x 



Dt V dt 



dx 
~dt 



Proof. Since v = f t {x,£ x ), ver(f) = v - hor(u), and hor(» = (x,£, x ,x, [t; x ,X(x)}), 
applying Definition 12 . 2 1 for the covariant derivative induced by a horizontal bundle, 
we conclude 



D(x,Q ( (d, \\ ( d& 



Dt " + \— \dt K ~''* B 'J J V dt ^ x,X \dt) 
as stated. □ 

The formula in Proposition l2.2l immediatelv yields the expression of the covariant 
derivative on the dual bundle, g Q , namely, 

D ( da x 

— {x, a x )= i x ,-jf- ad x(i) K 

This will be useful later when we take covariant derivatives of momenta. 

Proposition 2.3. Given an interpolation method, X : TX — > Xdi v (M), the map 
m x : T SDiff (M)/G e — >■ TX © 00 given by 

is an isomorphism of vector bundles. 
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m) 



Proof. We first must show that is well defined on the quotient space T SDiff (M)/G Q . 
Let ip G G G ; note that 

MK ° ^]) = KW©)), K VO o o v)- 1 - zOv(v>(©)))) 

= W0),^o^-%(0))) 

Additionally, it is easy to check that \&2 has the inverse \l/j 1 (a;, i, ^ x .) = [it o <p] 
where u = £ x + X(i). □ 

3. Lagrangian Reduction and the Equations of Motion 

The kinetic energy of a fluid flowing on M, denoted L : T SDiff (M) — > R, is 
given by 

(4) L{<p,<p) = % [ ||^(m)|| 2 rfvol( 

where p denotes the density of the fluid, assumed to be constant. It was shown 
in |Arn66j that L is SDiff (M)-invariant and that the resulting Euler-Poincare 
equations are precisely Euler's equations for an ideal, inviscid, homogeneous, in- 
compressible fluid 

du Vp 

— + y u u = 

ot p 

divw = 

for any oriented bondaryless Riemannian manifold M. In this section we shall 
take M = R d and will do a reduction by G Q C SDiff (M). Of course, the resulting 
equations of motion yield the same dynamics. However, the use of the interpolation 
method heavily influences how one writes the equations of motion. To see this, we 
must first study how variations of curves in SDiff (M) lead to variations of curves 
in TX © 00 under the map induced by an interpolation method, X. 



3.1. Covariant Variations. Let ip t be a curve in SDiff(M). Then a deformation 
is a surface embedding, <p\,t, such that <po,t = V 9 *- We desire to measure how much 
the variation 8tp t := changes the quantity (x, £ x ) := (x, <p o ip' 1 —Z(x)) G 

g Q where x(t) = v ? t(©) and (pt = -%r- To do this, we will invoke the covariant 
derivative induced by X. 

Additionally, through the Riemannian metric, (•, -) x , on X given by 



(v X} w x ) x := / (X(v x )(m),l(w x )(m)) M dvol(m), 
Jm 

we may use the resulting Levi-Civita connection to get a covariant derivative on 
the Whitney sum TX © q q . This direct sum of covariant derivatives will also be 
denoted 
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Using this covariant derivative, we may define the covariant variation of a curve 
(x,£ x )(t) G 00 with respect to a deformation (x, £ x )\t by 6 x (x,£ x ) := D ^'\ x ^ 

nA A=0 

However, we will primarily be concerned with variations induced by variations of 
curves in SDiff(M). The following propositions describe the form of such varia- 
tions. 

Proposition 3.1. Let tp t be a curve in SDiff(M). Set tp t = ±(<p t ) G T^ t SDiff(M), 
x(t) = (f t (Q),x(t) = (fi t (Q) = f . Finally, set 

-l 



£ x = ipoip -l(x), 

so that we have a curve (x,£ x )(t) G g Q . Then, provided a vertical deformation of 
<f t given by p t ,\ = Pt ° i>t,x f or a deformation i/j t ,x G G Q where ip t ,o ^ the identity, 
the covariant variation of (x,£ x )(t) is given by 



where i] x = (p*r) and rj := 



A=0 



Proof. We decompose 5£ x : = ^| A=0 (£ x ) into three parts: 

dip i dip^ 1 dl(x) 



Ti T 2 T 3 

As X(x) does not depend on A we may set T 3 equal to 0. We can rewrite 7\ as 



•9 



/ • \ -i d 2 ip _! 



= ^(Tfori)oip 1 

= Tip O 1] O p~ l + T(/9 O f] O p 
= T(pOT] O ip' 1 + ^7), 

while T 2 may be written as 

T 2 = T(po = -Tip o v?- 1 o 5p o <^" 1 

CM 

= —Tip O o p o r] o p 

= —Tip O 7] O p^ 1 

Thus, we find 5^ = p>J). Since <5x = 0, we see that 

$ X (x,£ x ) = (x,8£ x + [I(6x),£ x ]) = (x,(p*rj). 



-i 



10 HENRY O. JACOBS, TUDOR S. RATIU, MATHIEU DESBRUN 

Additionally, note that 
We calculate 



= -[£ + X(x),r) x ] + 
Therefore v?*?7 = -j- + [£ + X(x),rj x ] so that 

5 x (x,e.) = (x, ^ + [e e +x(i),%]) = + [(x.6,), 

as required. □ 

While vertical variations are now clear, we must also consider how (x, £ x ) = 
(x,(p o tp^ 1 — X(x)) varies in response to variations of x £ X. Given a curve 
(ft £ SDiff(M) we may take a deformation of the curve x t = <pt{&) £ X given by 

x t) \. If we consider the sequence of vector fields X {^-^f^ parameterized by A we 

may integrate them over A do get a diffeomorphism. This defines the horizontal 
deformation of ip t by 



(5) <Pt,x = / X 



x ' dx 







f.A 



<9A 



d\' o yj t . 



Proposition 3.2. Let tp t be a curve in SDiff(M) and let x(t) = <p t (Q). Let x tj \ be 
a deformation of x{t) and let tp t ,\ be the resulting horizontal deformation given by 
equation fl5]). Then the covariant variation of (x,£ x ) where £ x := if o ip~ Y — X{x) 
is given by 

5 x (x,£ x ) = B(x,5x), 

where 

(6) B(x,5x) = (x,dX(x,5x) — [X(x),X(8x)]). 

Proof. Upon invoking equation ([Tj) and the fact that the fiber derivative, ||, is 
identical to X (because X is linear on fibers of TX) we find that 



o 

A=0 



ip l -ipoip 1 o5ipoip 1 - (^-(x),5x^ -X (J^.( 5x )J ■ 
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Using equality of mixed partials we find that 

dip 



0i P 1 = -n(S<p)°<p 1 

A=0 0,1 



— (8x),x\ +1 ( £^(<fo)J -X{SX) OLpOLf- 1 



Upon substitution into the last line of the previous calculation we find 

5£* = ((^(8x),x^ - (j^;{x), 6 x yj + X{5x) o (0 o ^ + (0 o ^) o X(5x) 

= dl(x, 5x) + [I(8x), tp o (p^ 1 ] 
= dl(x, Sx) + [l(8x),£ x + Z{x)\. 

Therefore, the covariant variation is 

8 T (x,£ x ) = (x,8£ x - [l(5x),£ x ]) = dl(x,8x) - [l(x),l(5x)] := B(x,8x) 

as stated. □ 

We call B the reduced curvature tensor. It measures the non-integrability of the 
horizontal distribution induced by the interpolation method. This becomes even 
more clear in the following proposition. 

Proposition 3.3. The reduced curvature tensor fl6]) may equivalently be written 
as 

B{x,8x) =l([l(x),l(8x)](x)) - [X{x),X{8x)}. 

The proof is built on the standard method used for the curvature tensor asso- 
ciated to a principal connection. In order to preserve the flow of the presentation 
we have included the proof in Appendix [A] 

As a result of Propositions 13.11 and 13.21 we see that the most general covariant 
variations of (x,^ x ) G 0© that are needed in the variational principle, are of the 
form 

(7) 8 x (x,Q = + [(*,£*), (x, Vx )} + B(x,8x) 

for some curve (x,r] x ) G g© and a variation 8x(t) of the curve x(t). 

3.2. Lagrange-Poincare Reduction. In this section, we state the Lagrange- 
Poincare reduction theorem (see [CMR01J) in terms of interpolation methods 
rather than principal connections. The resulting equations of motion induce the 
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same dynamics, of course. Note that the reduced Lagrangian, [L], written on 
TX © g via the isomorphism is given by 

(8) [L](x,x,Q = ?- [ ||£, + X(x)|| 2 dvol(m). 

Theorem 3.1. Let L : TSDiff(M) — > R be a G Q invariant Lagrangian such as 
the kinetic energy Lagrangian given in equation (p0). Let I : TX — > 3£div(^0 be an 
interpolation method, and let [L] : TX © $j — > M fre i/ie reduced Lagrangian. Let 
(ft G SDiff(M) 6e a curve and sei x(t) = <ft{0)> £x — <Px ° V 9 " 1 — 1(x). Then the 
following are equivalent. 

(1) The curve <p t is critical for the action 

S = L((pt,(p t )dt 



with respect to variations 5(p t with fixed end points. 
(2) The curve tp t satisfies the Euler- Lagrange equations 



Dt \d(p J dip 

with respect to an arbitrary covariant derivative and connection on T SDiff (M) . 

(3) The curve (x,x,£ x )(t) is critical for the reduced action 

[S] = / [L](x,x,^ x )dt 
Jo 

with respect to arbitrary variations 5x(t), with fixed end points, and covari- 
ant variations of (x, £ x ) of the form 

S x (x, £,) = - [(*, V x), (x, Q) + B(x, 5x) 

for some curve (x,r] x ) G . 

(4) The curve (x, x, £ x ) satisfies the Lagrange- Poincare equations 

urn -«•«<» 

where B^{y X) w x ) = (^^, B(v x ,w x ] 

Proof. The equivalence of (1) and (2) is an intrinisic formulation of the standard 
derivation of the Euler-Lagrange equations such as done in |AM00j . We have 
chosen to write the intrinsic formulation more out of neccesity than interest, as 
we are working on a space with non-trivial coordinate charts. In this case, the 
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"equivalence of mixed partials" comes from our definition of the covarient deriva- 
tive induced by a generalized connection. Specifically, let (p t>e be an embedding of 
a surface into SDiff(M) (i.e., a deformation of a curve). Then we observe that 



(9) 

where <p t 



dt 



De 

and Sift 



v l ver 



dipt 



d 2 cp 







D6<p t 
Dt 



. Using this, one is able to prove the equiva- 

=o 

lence of (1) and (2) in the standard way using integration by parts. 

We now prove the equivalence of (1) and (3). If the time integral of L is ex- 
tremized along (ip,ip)t then [L] must be extremized along ty x ([<p t ,<p t ]). However, 
by proposition 13.21 and 13. 1[ we know that it must be extremized with respect to 
variations of the form given in equation (J7J). 

Finally, we prove the equivalence of (3) and (4). Assume [S] is extremized with 
respect to the variations given in equation (J7J). Then we find, with equation (QQ), 
that 



5[S) 



d_ 

dX 

t 



A=0 JO 



[L](x x ,x x ,t; x ,\)dt 



(d[L],8(x,x,£ x )) dt 



/d[L) 



\d£ x 



8% 



\ dx 



dt. 



/d[L] Dx\ 

Using the definition of the covariant derivative and equality of mixed partials we 
get 

Dx ( ( d 2 x 



De 



ver 



V dtde J 



D5x 



Dt 



so that the last calculation becomes 



5[S] 



/d[L] D(x, Vx ) 



+ 



\ d£ x ' Dt 
/ d[L] D5x 
\di'~Dt 

D_(d\L[ 
Dt V d£ x 

D_(m 

Dt I dx 



+ [(x, r] x ), (x, £ x )} + B(x, Sx) ) dt 



+ 



\ dx 
— adt 

dx 



5xj dt 



dt 



i±B^5x 



dt. 



The theorem follows from the arbitrariness of (x, r] x ) and Sx on the interior of the 
interval. Additionally, the above sequence of calculations is reversible. □ 

At this point, one may be inspired to come up with particle methods by trying to 
better understand the horizontal equation (the equation for the dynamics on X). 
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Alternatively, others may be very skeptical of this idea since the vertical equation 
is nearly everything. In fact if we unpack the terms of the vertical equation, we 
find that = (£ x +X(x), •) L2 = u b , so that the left hand side is 

while the right hand side of the vertical equation is — adL^s (u ) . Bringing both 

terms to one side we find %■ + ad* u b = 0, which is the inviscid fluid equation 
|AK92j . except for the fact that the vertical equations (strictly speaking) only 
address the domain M\x. In essence, the horizontal equations only state that the 
particles move in such a way that u can be extended smoothly by "filling the hole" . 
Despite this sobering observation, we know that computational scientists simulate 
fluids and successfully use interpolation methods all the time. After studying 
what happens when we reduce by a class of subgroups of G© in §HJ we will try to 
understand how the horizontal equation can potentially inspire particle methods 
in 83 

4. Higher Order Isotropy Groups 

In the previous section we reduced our system by the Lie group G . In this 
section, we consider the Lie group 

G$ := {i/j G G© | T*ij) is the identity on T*M}. 

In local coordinates, elements of Gq are diffeomorphisms such that the Taylor 
expansion around each point of has a 1 for the first coefficient, a for the 
coefficients 2 through k, and is arbitrary for the remaining coefficients. To make 
this more precise we will include a short discussion on jet bundles — in particular, 
jets of elements in SDiff(M). 

4.1. Jet Bundles of the Special Diffeomorphism Group. Consider the fol- 
lowing equivalence relation: ipx, (p 2 £ SDiff (M) are equivalent to kth order at if 
they have the same Taylor expansion at 0. We denote the set of equivalence classes 
for this relation by j£ (SDiff (M)). Elements of J"* (SDiff (M)) are called fc-jets of 
SDiff (M) sourced at 0. For k = this simply means that (pi((D) = ^2(0) and for 
k — 1 this means T & ipi = T & ip 2 . Given a G SDiff (M), we denote its fc-jet by jj^ip. 
Additionally, there is a natural projection vrf : j|(SDiff (M)) ->■ J"|(SDiff(M)) for 
I < k. In particular, the projection 7Tq sends fc-jets to iV-tuples of points in M at 
the base of these jets via the map j (<£>) | — > V 9 !©) e X^°\ 

In addition, "jet" is a functor and it commutes with the "tangent functor T". 
This means 

TJ"^(SDiff(M)) = Jq(T SDiff (M)). 

Given a vector-field u G Xdi v (^) and an o^ -* G X, we may consider the fc-jet 
j k (o)( u ) by viewing a as a tangent vector to SDiff (M) at the identity. Given 
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any ip E SDiff(M) such that <£>(©) = we see that j%(u o tp) is an element of 
TJq (SDiff(M)). Moreover, we may consider the entire equivalence class, x = j^ip, 
as a set of maps which can be composed with the set of maps comprising j fe (0) (u) 
to arrive at j k (0) (u)ox = j^(u o ip). The notation used on the left hand side of this 
equation will be useful when comparing ideas in this section to those mentioned 
in the previous sections. 

Finally, we can also consider the set J~q(G q ) which consists of equivalence classes 
restricted to elements of G & . It should be clear that Jq(G & ) is a (finite dimen- 
sional) Lie group with the group multiplication j%(ipx) ■ i%(ip2) '■= j'o^i °i ) 2)- The 
Lie group Jq(G q ) acts on j|(SDiff(M)) through the right action 

for each ip E SDiff (M) and ip E G & . We will use this action in the context of 
Noether's theorem in §5] to find conserved momenta for a class of particle methods. 

4.2. Quotients and k th order Interpolation Methods. In this section we 
describe the quotient SDiff {M)/d§. 

Proposition 4.1. The quotient X^ := SDiff (M)/G ( q ) is given by k-jets ofSDiS(M) 
sourced at 0, i.e., 

= Jl (SDiff (M)). 

Proof. Let ip E SDiff (M). Then we see that j^ip = j^ip o ip) for all ip E G k & . 
So, for each equivalence class [<p] there exists a unique fc-jet. Conversely, a k- 
jet x E (SDiff (M)) can be associated to the equivalence class of some <p E 
SDiff (M), which is only unique up to an element of G@, by the very definition of 
J| (SDiff (M)). □ 

Example 4.1. Upon choosing a frame {fi, ■ ■ ■ , f n } with unit volume over the 
points ©i, . . . , Oat, we may identify X^ as the set of iV-tuples of unit volume 
frames with non-overlapping base points. If M = M d we can view unit volume 
frames as elements of SL(d) and we can view as the trivial fiber-bundle over 
given by the set X^ ^ x SL(c?) 7V . Through this observation It should be clear 
that X^ consists of particles which carry eccentricity and orientation. For k > 1 
we may interpret fiber bundle over X^ where the fibers consist of 

infinitesimal perturbations of orientation and shape. Even less formally, X^ can 
be visualized as iV-tuples of non-overlapping particles with orientation, shape, and 
an ability to "jiggle" like rubber to k th order. 

Let us return to the general case and follow the procedure of Lagrange- Poincare 
reduction used in the last section but for a kth order interpolation method, formally 
defined below. 

Definition 4.1. A k th order interpolation method is a map, T : TX^ k ' — > Xdi v (M), 
such that j k {0) (1(A)) o x — (x, A) where (x, A) E TX^, and x^ = Hq (x). 
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In other words, a k order interpolation method takes infinitesimally varying 
fc-jets (given by x) to vector fields which match these infinitesimally variations 
fc-jets at the points of the particles (given by x^). Even more concretely, k th 
order interpolation methods approximate vector-fields to k th order, given k th order 
constraints at a finite set of points. 

By using a k th order interpolation method, we can follow the same sequence we 
used before by defining the isomorphism ^ x : TSDiff(M)/G'5 ) -»■ TX<® © g$ 
given by 

*z(K}) = © K ° v -1 - X M) ■ 

From here the symbolic manipulations are entirely the same as in the previous 
sections. The reduced Lagrangian 

[L](x,x,Q = - I \\X(x) + £, x \\ 2 dvol(m). 

1 J M 

is now defined on the Whitney sum, TX^®q§. The La grange- Poincare equations 
for the reduced Lagrangian system are described on this Whitney sum as well. The 
horizontal equation, 

Dt \ dx J dx x M 

on TX( k \ describes the motion of particles with extra structure (corresponding to 
fc-jets of diffeomorphisms) , and the vertical equation, 

on fibers of the associated bundle, is formulated in terms of vector fields which 
vanish to k th order at the particle locations. Additionally, the vector field I(x) 
approximates the spatial velocity u = (p o ip~ x to k th order around the particles. 

Remark 4.1. There is an extra symmetry in the horizontal equations. In particular, 
the group Jq{G q ) is a remaining symmetry left behind by G@. If M = M. d and 
k = 1 we can identify Jq{G q ) with SL(<i) Ar by choosing unit volume frames above 
©i, •• .,0jv- 

If the Lagrangian does not depend on the orientation and shape of the particles, 
we get a right SL(<i) Ar symmetry which has yet to be "quotiented away". This 
extra symmetry is attached to the particles, and results in conserved momenta 
as a consequence of Noether's theorem. We will discuss this more in §5.21 in the 
context of a particle method. 

5. Applications to Particle Methods 

We hope that a computational scientist or engineer interested in particle meth- 
ods for fluids finds the Lagrange-Poincare equations thought provoking; we have 
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deliberately presented our approach using the "interpolation" point of view to en- 
courage such thoughts from practitioners. In this section we describe two methods 
by which one can construct particle methods. This is not to say that these are 
the only methods one can do! Study along these lines is wide open for further 
exploration and implementation. 

5.1. A Variational Method. The horizontal equations look very similar to the 
Euler-Lagrange equations except for the fact that there is a curvature force on the 
right hand side and the particle momenta, depend on the vertical components 
(x,£ x ) G 0Q. It is tempting to find a way of ignoring £ x . One way to do this is 
to add a fictitious non-holonomic constraint which restricts the spatial velocity 
u = ip o ip" 1 to be in the range of X. That is to say, we introduce the constraint 
distribution V x C T SDiff(M) by defining its fibers; the fiber over if G SDiff(M) 
is 

V x := {X(x) o ip | x G T x X<- k \x = j*<p}. 

Then we restrict L to T> x . By construction, the non-holonomic force which keeps 
Lp in V x is such that the vertical component £ x = Lpo^p" 1 — T^ k \x) vanishes, where 
x = <£>(©) • Thus the dynamics on will determine the dynamics completely. 

By construction, the constraint distribution V 1 is Gq ^-invariant. This means 
that for each v v G V 1 , it follows that v v o G T> x for all ip G Gq\ We denote 
the set of these orbits by [D 1 ]. Due to the injectivity of X, there exists a vector 
bundle isomorphism $x : TX^ — > [D x ] which sends x to X(x) o x where we view 
x = T X (k) (x) as an equivalence class of diffeomorphism, so that it makes sense to 

(k) 

compose it with another map. Given a G Q -invariant Lagrangian, it should be 
clear that X induces a Lagrangian on TX^ given by L x = [L]\ vX o& x where [L] is 
the Gq ^-reduced Lagrangian. We may use this construction to define a Lagrangian 
system on as in the following theorem. 

Theorem 5.1. Let L : TSDiff(M) R be a Gg° -invariant Lagrangian and let 
ip t be a curve in SDiff(M). Let T> x be a velocity constraint distribution induced by 
a kth order interpolation method, X. Finally, let ip t be a curve in SDiff(M) and 
set x(t) = j^ift- If [L] is the reduced Lagrangian and 

(10) L x := L\ vx o $ x : TX™ -)• K, 

then the following statements are equivalent: 

(1) (fit satisfies the constrained Euler-Lagrange equations with respect to the 
velocity constraint distribution V x . 

(2) The unreduced action S = J* * L(ip, ip)dt is extremized with respect to fixed 
end-point variations of the curve ip t in restricted to the distribution V x . 

(3) The action Sx = Jq Lx(x, x)dt is extremized with respect to arbitrary vari- 
ations of the curve x(t) with fixed endpoints. 
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(4) The curve x{t) satisfies the Euler- Lagrange equations on with respect 
to Lx- 

Proof. Assume S is extremized with respect to variations in T> x . In other words, 
S J L(ip, ip)dt = with respect to variations Sip = X(Sx) o ip for some variation Sx 
of the curve x(t). This works, because each fiber T> x is isomorphic to the tangent 

fiber T x X^ k \ Since both L and T> x are G@ ^-invariant we find 

S J L{ip,ip)dt = S J [L]([ip,ip])dt, 

where the variations on the right hand side are taken in [D x ] by using [Sip]. 
These variations are given explicitly by $(<5x). Furthermore, the isomorphism 
$2 : TX (fc) -> [V 1 ], allows us to verify the equations L(ip,ip) = [L]($i(x)). How- 
ever, the right hand side of this expression is Lx- Therefore, extremizing Sx with 
respect to arbitrary variations of x(t) is equivalent to extremizing the restriction 
of S to T> x with respect to variations in T> x . The equivalence with the constrained 
and unconstrained Euler-Lagrange equations is a well known result. □ 

Theorem 15.11 suggests that we can estimate the exact evolution of an ideal fluid 
by solving the Euler-Lagrange equations for the Lagrangian Lx — [L] o $ on the 
space and then use the interpolation, X, to approximate the spatial velocity 
field. This is a method for estimating the flow of an inviscid fluid using only particle 
positions and velocities, i.e., a particle method. The error of such a scheme may 
be measured by the magnitude of the Lagrangian parameter used in constraining 
the dynamics to T> x . 

Remark 5.2. If M = M. d , the constraint force to V* x is the (co)vector field 1.[x) ■ 
V(X(x)). Noticing that this is computed solely in terms of particle motion, suggests 
that we can take a norm of this quantity and use Grownwall's inequality to find 
an error bound for our particle method. Additionally, this quantity is generally 
a function solely of x and, therefore, can be used as a stopping criterion in these 
integrators. 

5.2. Kelvin's Circulation Theorem for Particles. Recall the remaining sym- 
metry of the horizontal equation mentioned in Remark 14.11 If the k th order in- 
terpolation method is compatible with this symmetry, then the numerical method 
suggested in Theorem 15.11 will also possess this symmetry. This is a valuable 
property for users who desire to preserve this qualitative property of the system. 
Additionally, it allows one to construct integrators that conserve momenta (via a 
discrete time Noether theorem); see |MW01] or |HLW02] for details. We formalize 
this statement in the following proposition. 

Proposition 5.1. Given the constraint distribution, V> x , induced by a k th order 
interpolation method, I, the quantity 

J c (x) = (x,x o£) x 
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is conserved by the Euler Lagrange equations for the Lagrangian, Lx, defined by 
equation ([TO]) for each £ in the Lie-algebra of J~ k (G & ) . 

Proof. As already mentioned, Lx is right invariant with respect to the J k {G Q )- 
action. The action of g G J k {G®) on an x G X k is given by composition of jets, 
x i — y x o g. Thus, the infinitesimal generator of £ is given by the vector field 

Therefore, Noether's theorem implies the desired result upon noting that the 

ILx 

dx 



momentum, is given by the one-form (x, •) x . □ 



/ " * * : < U 

M * * T f 

Figure 2. A particle-centric version of Kelvin's circulation theorem 
corresponds to conservation of circulation along infinitesimal curves 
circling the particles. For k — 1 this implies that the spin and the 
rate of stretching are conserved. 



Proposition 15.11 may be viewed as a particle-centric version of Kelvin's circu- 
lation theorem. In particular, it is known that the conserved quantity for the 
particle relabeling symmetry of an ideal fluid is given by the circulation. This 
implies that the integral of the dot product of the fluid velocity with an arbi- 
trary vector field along an advected closed curve is conserved by the inviscid fluid 
equations |AK92} Chapter 1]. Proposition 15.11 tells us that curves which are in- 
finitesimally close to the particles (see Figure H]) conserve the same quantity, even 
for the non-holonomically constrained version of the inviscid fluid equations de- 
scribed in Theorem 15.11 If A; = 1, this implies that the flow will conserve the 
spin and the rate of deformation of each of the particles. For k > 1 we can inter- 
pret the conserved quantities as manifestations of Kelvin's circulation for k th order 
perturbations of these infinitesimal curves circling the particles. Therefore, when 
modeling the fluid using the finite dimensional Lagrangian system on with 
Lagrangian Lx, one can better control the conservation of momenta and Kelvin's 
circulation theorem. 

Remark 5.3. Variational particle methods for the EPDiff equation have been imple- 
mented and analyzed in |CDTMl2] using the reconstruction mapping mentioned 



20 



HENRY O. JACOBS, TUDOR S. RATIU, MATHIEU DESBRUN 



in Remark 12.11 The convergence of this method was used to prove global exis- 
tence and uniqueness in |CLP12j . Therefore, these ideas are not totally without 
precedent. However, the cited references are restricted to the case where k = 0. 

5.3. Spectral Methods and Hybrid Spectral-Particle Methods. Building 
upon the ideas of the previous section, we can consider a method which is in- 
dependent of the particle locations. To do this, choose N linearly independent 
vector fields, Ui,...,Un G Xdiv(M). Then construct the interpolation method 
X : TX^ Xdiv(M) given by 



where the coefficients q are described implicitly by the constraints of definition 
14.11 (it is a linear algebraic inverse problem). Using this interpolation method pro- 
duces a spectral method since solving the Euler-Lagrange equations is equivalent 
to extremization over the constant vector space V = span(iti,--- ,Un) at each 
time step. This maybe counter intuitive, since spectral methods do not depend 
on the location of particles. One clear way to see the proposed method is not 
dependent on the locations of the particles is to note that given two different par- 
ticle configurations x and y in X^°\ the velocities x = CiUi(x) and y = CiU^y) are 
both interpolated to the same velocity field, i.e., I{x) = I(y) = c(iii. The method 
obtained by solving the Euler-Lagrange equations for the Lagrangian Lx defined 
in equation ( Tl0|) . is equivalent to choosing the coefficients c$ at each time in a 
manner guaranteeing that the action is extremized. In discrete time a variational 
integrator would generate the constants q at each time step in such a that they 
do not depend on the particle positions; but only on the previous group of q. 

Remark 5.4. One should probably choose the basis {ui, . . . ,un} in such a way 
that as the number of particles goes to infinity, the basis captures a function space 
of decreasing regularity, e.g., using a Fourier basis on M = M. d . 

Judicious choices of bases lead to nice error bounds. One does not need to 
choose a full basis, however. One could choose vector fields for k < 

N to construct an interpolation method which is partially spectral and partially 
dependent on particle locations. This can be considered a "multiscale" model in 
the sense that the spectral basis keeps track of the "large scale dynamics" while 
the remaining degrees of freedom are described explicitly by the motion of particles 
and can be thought to model the "fine scale dynamics" . 

5.4. Vortex Methods. It was shown in iOSOll that Chorin's vortex blob algo- 
rithm (originally developed for simulating inviscid fluids) provides an exact solu- 
tion to the equations of motion for am ideal, inviscid, homogeneous, incompressible 
second grade fluid. More specifically, Chorin's vortex blob algorithm yields solu- 
tions to the Lagrangian system on SDiff(M) with the reduced Lagrangian given 
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in terms of the spatial velocity by 

[L] a (u) := - l (|H| 2 + a. ■ trace (Vw T Vw)) d 2 x, 

for some a > 0. By construction, the vortex blob method conserves circulation 
by advecting guassian blobs of vorticity. In this section we will describe in words 
how the vortex blob method may be viewed as a particular case of the algorithm 
described in this paper. Consider a I s * order interpolation method, X, such that 
the spin of each particle is mapped by X to the spatial velocity field of a vortex 
blob, and the translational component of the particles is mapped by X to a spatial 
velocity field with vorticity above each particle. Then, by Proposition 15.11 if we 
initialize the system with a spatial velocity field consisting solely of vortex blobs 
located above the particles, it follows that X(x) remains a sum of vortex blobs for 
all time. Additionally, the particles will be advected by these vortex blobs and the 
vorticities will extremize the action of [L] a . As mentioned at the beginning of this 
section, this is precisely the evolution of the vortex blob algorithm. 

Additionally, there are variants of the vortex blob method which allow for de- 
formations of vortex blobs, as would be the case of an algorithm induced by a 
k th order interpolation method for k > 1. For example, |UWB12] uses Hermite 



functions to generate a basis of functions above each vortex blob to do just this. It 
is not immediately clear whether this method satisfies a variational principle and 
therefore it is not known whether such a method matches the framework described 
in this paper. Nonetheless, appealing to higher-order deformations of vortices did 
yield higher accuracy in numerical experiments for the case k = 2. Another in- 
teresting method potentially linked to our approach, is the adaptive smoothed 
particle hydrodynamics proposed in [SMV096J. 



6. Extensions to other fluids 

Given the structure presented so far, the process of extending these ideas to 
other types of ideal fluids can be summarized with little effort. In this section, we 
will provide a sketch of how to extend the constructions of the previous section to 
complex fluids and Euler-a models. 

6.1. Complex Fluids. Consider the set J-(M, O) of smooth maps from a mani- 
fold M to a Lie group O. The Lie group SDiff (M) acts on J-(M, O) on the right 
by the assignment / 6 J*(M, O) ^ / o cp e 7{M.O) for each Lp e SDiff (M). 
This right action defines the semi-direct product Lie group SDiff (M) (s)X"(M, O), 
where the group multiplication is (ipi, fi) ■ (y?2, $2) := (^1 °V ? 2, (/i ^) • fi)- A uni- 
fying framework for ideal complex fluids was described in [GBR09J by performing 
Euler-Poincare reduction on a Lagrangian system with the configuration manifold 
SDiff (M) (§)J-(M, O). The number of scenarios captured by this framework is 
truly remarkable (magnetohydrodynamics, spin-glass fluids, liquid crystals, fluids 



22 



HENRY O. JACOBS, TUDOR S. RATIU, MATHIEU DESBRUN 



with microstructure, Yang- Mills fluids, etc.). It should be clear that one may re- 
duce such a system by Gq ' as well. This could be particularly fruitful for the case 
of fluids with microstructure, where the local orientation and shape of the fluid 
particles shows up in the Lagrangian. For example, in the case of micromorphic 
nematic liquid crystals, O = SO (3). We see that for M = 1R 3 is the triv- 
ial bundle (R 3 ) N x SL(3) JV . Since SO(3) is a subgroup of SL(3) it follows that a 
method on X^- 1 ' could be capable of modeling liquid crystals. 

6.2. 2nd Grade Fluids and other generalizations. A fairly modest general- 
ization to the preceding setup is to add gradient terms to the Lagrangian and 
consider second grade fluids. For example, the kinetic energy on M = M. 3 given by 



where u = (p o ip^ 1 and a > 0, clearly exhibits particle relabeling symmetry and 
leads to the Euler-a model in the incompressible case. We may even consider the 
Lagrangian 



where I is a positive definite differential operator on X(M) (again, the Helmoltz 
operator I — a ■ V is a good example). The momentum is given by G * u where 
G is the Green's function for I. In this case, we may even consider compressible 
fluids if we use interpolation methods which map to the full set of smooth vector 
fields on M (see jl 1SS09. Chapter 11] for a good overview). If G is non-singular, 
the interpolation method mentioned in remark 1231 becomes a particularly natural 
choice. 



Since the publication of [Arn66j, it has been well known that Euler's equations of 
motion for an ideal, inviscid, incompressible, homogeneous fluid may be regarded as 
Euler-Poincare equations on the Lie algebra of divergence free vector fields tangent 
to the boundary. Therefore, it has been clear that one could reduce by various 
subgroups of the diffeomorphism group. In this paper we have accomplished one 
version of this reduction, using the notion of an interpolation method. We chose 
to use interpolation methods instead of principal connections (as in [CMROlj ) 
to clarify the notion of estimating the velocity field of the fluid with particles. 
Specifically, if an interpolation method is chosen, one may estimate the spatial 
velocity field of the fluid and even estimate the evolution of the system over short 
times by integrating a non-holonomically constrained version of the equations. 
We also discussed the reduction process for higher-order isotropy groups. The 
horizontal equations for higher-order isotropy reduction contain extra symmetry. 
Using a particle method which respects this symmetry produced a discrete version 
of Kelvin's circulation theorem. 





7. Conclusion 
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Future work includes: 

• Implementing an instance of the particle method described in this paper 
for k > 1 and comparing performance with the vortex blob method. 

• Extending these constructions to other types of fluid, particularly complex 
fluids. 

• Appending a finite dimensional model of the vertical Lagrange-Poincare 
equations would yield a finite dimensional model of a fluid on a Lie groupoid. 
This could generalize the integrator described in |GMP + 11] which models 
fluid on a finite dimensional Lie group. 

Appendix A. Interpolation Methods and Principal Connections 
In this appendix our goal is to prove that 

B(x,6x) =X([X(x),X(6x)](x)) - [X(x),X(5x)\ 

where B is the reduced curvature tensor defined in equation In the process, we 
will also clarify the "equivalence" of interpolation methods with principal connec- 
tions on the right principal bundle it : SDiff (M) — y X. First, we use the definition 
of a (right) principal connection (as in [MMR90J, but adapted for the case of right 
invar iance). 

Definition A.l. Let it : Q — )■ S be a right principal bundle with structure group 
G. Let q be the Lie algebra of G. A right equivariant principal connection on Q 
is a Q-valued one-form A G Q X (Q, q), satisfying the following properties: 

(1) For each £ G g, we have A(£q) = £, where £q e %(Q), defined by ^(q) : = 
~^] t - Q ' ex P(^0 f or an V q & Q, is the infinitesimal generator o/£. 

(2) For each g G G and v G TQ we have A(g ■ v) = Ad~ 1 (A(v)). 

In the context of the principal bundle SDiff (M) over X, the notion of a right 
equivariant principal connection is identical to the notion of an interpolation 
method. This is shown in the following lemma. 

Lemma A.l. Given an interpolation method X : TX <^-> Xdi v {M) the map 
(11) A(v v ) = T^ 1 o Vip - V *X{v v {Q)) 

is a right equivariant principal connection on the G & -principal bundle ttx '■ SDiff (M) — y 
X. 

Proof. We must prove three things. 

(1) Show that A maps to the Lie algebra of G & . 

(2) Show that A satisfies the first property of a principal connection. 

(3) Show that A satisfies the second property (i.e., right equivariance). 
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To show that A maps to the Lie algebra of G & is equivalent to proving that 
the values of A are divergence free vector fields vanishing at all points of © C 
M. Clearly, Tip' 1 o v v is a divergence free vector-field on M, as it is merely 
the left trivialization of v v , and ip*I(v v (o)) is a divergence free vector field by 
construction. Thus, the difference of the two, namely Aiv^) = Tip" 1 o v v — 
ip*(Z(v v (Q))), is a divergence free vector field as well. Additionally, we have 
A(v v> )((D) = Tip" 1 (v v ((D)) — Tip" 1 (v v (&) = 0, where we have used the property 
I(v x )(x) = v x for all v x G T X X. This proves (1). 

To prove (2), let £ G g . The infinitesimal generator of £ is the map ip h-> T<po£. 

It is now easy to check that AiTip o £) = £ for all £ G $j . 

Finally, to prove (3), note that for any ip G G Q and v v G TSDiff(M) we have 

A( Vip o V) = T(cp o V)" 1 o ( Vv oiJ>)-((po ip)*Z((v v o ^)(©)) 
= ^*(T^- 1 o^-^X(^(0))) 
= Ad^ 1 (A(v v )) 

since the adjoint action is given by push-forward. □ 

Conversely, an arbitrary right equivariant principal connection, A G Q 1 (SDiff(M), g Q ), 
naturally induces an interpolation method. This is done as follows. Recall that 
the kernel of A is a horizontal subbundle of TSDiff(M) which naturally defines 
the horizontal lift operator : TX © SDiff(M) -»■ TSDiff(M). An interpolation 
method, X, is obtained by composing this horizontal lift with the right trivializa- 
tion of TSDiff(M), i.e., " ' 

= X(x) o ip. 

The above equation defines X uniquely. Moreover, the principal connection ob- 
tained from X in equation pip is precisely A. Thus the concepts of principal 
connections and an interpolation methods are effectively equivalent for the princi- 
pal bundle vr : SDiff(M) -+ X. 

It is shown in [MMR90J that the reduced curvature tensor ([H]) may also be given 
by the expression 

B(x,6x) = v.A([hl(x),hl{6x)]), 

where the choice of ip is arbitrary up to the constraint y?(0) = x = tx(x). This 
version of the reduced curvature tensor allows us to prove Proposition 13.31 in view 
of the following lemma. 

Lemma A. 2. Let [,]m be the Lie bracket of vector fields on M and [,]sDiS(M) the 
Lie bracket of vector fields on SDiff(M). Let p triv : TSDiff(M) — > j£di V (M) be the 
right trivializing morphism 5ip i— >■ 5ip o ip" 1 . Then, for X, Y G X(SDiff (M)) ; we 
have 

[Ptriv ° A, p triv oY] M = p triv o [A, ^] SDiff ( M ) • 
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Proof. Using the dynamic definition of Lie derivative and evaluating at a given 
diffeomorphism <p £ SDiff(M) we have 

[X, F] S Diff(A/)(v 2 ) = d t d s (p t>s - d s d t ip t)S 

where d s = £\ S=Q , d t = f t \ t=Q and ip s>t is such that y? ,o = V, d t ip t ,o = X{ip) and 
d s (p , s = Y(cp). Applying p triv , we find 

Ptri Y {[X, F] S Diff(M)(^)) = d t d s ((p t>s o ip- 1 ) - d s d t ((p ttS ° y' 1 ) 

= [p tTiv (X((p)),p tTiv (Y((p))] M 
thus proving the statement. □ 
The take-away message from this lemma is that for each ip £ SDiff(M) 
[X, Y} SDiS(M) (ip) = [X{(p) o ip- 1 , Y(jp) o (p~ x \ M o ip 

for X, Y £ X(SDiff(M)). The right-hand side only depends on the values of X 
and Y at cp £ SDiff(M) and provides a bracket on the vector space T v SDiff (M), 
as opposed to the set of vector fields j£(SDiff(M)). With Lemma [A. 2| we are now 
prepared to prove Proposition 13.31 

Proof. Using Lemma IA.2I we may free ourselves of the burden of extending vectors 
in SDiff (M) to vector fields. This allows us to make the calculation of the reduced 
curvature tensor as follows. 

B(x,5x) = <p*A([hl(x),hl(8x)\) 

= p^Tip' 1 o [Z(x) o (p,X(5x) o(p]— l([l(x),l(8x)](x))) 

= [l(x),l(5x)]-l([l(x),l(5x)](x)) 

which proves the proposition. □ 
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